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Abstract - A spin- 2 Kitaev sublattice interacting with a subsystem of spinless fermions is studied 
on a honeycomb lattice when the fermion band is half filled. The model Hamiltonian describes a 
topological Kondo lattice with the Kitaev interaction, it is solved exactly by reduction to free Ma- 
jorana fermions in a static Z 2 gauge field. An yet unsolved problem of a hybridization of fermions 
and local moments in the Kondo lattice at low temperatures is solved in the framework of the 
model proposed. The Kondo hybridization gap is opened and the system is fixed in insulator 
and spin insulator states, due to a spin-fermion nature of the gap. We will show that the hy¬ 
bridization between local moments and itinerant fermions should be understand as hybridization 
between corresponding Majorana fermions of the spin and charge sectors. The RKKI interaction 
between local moments is not realized in the model, a system demonstrates a ‘quasi Kondo’ sce¬ 
nario of behavior with realization chiral gapless edge states in topological nontrivial phases. The 
ground-state phase diagram of the interacting subsystems calculated in the parameter space is 
rich. 


Introduction. — The Kondo lattice problem remains 
m unsolved problem of strong correlated systems, with 
nteracting delocalized and localized electrons m The 
£ondo lattice Hamiltonian describes conduction electrons 
nteracting with local moments arranged regularly. The 
<^ondo lattice approximation is used for description of 
■are-earth and transition compounds. The scenarios of 
lybridization of electrons and local moments, opening of 
;pin and charge gaps at low temperatures, a formation of 
arge volume of the Fermi surface in the Kondo lattice are 
;till unclear. In the framework of ID model of the Kondo 
attice Tsvelik has shown that the insulating state forms 
lot due to a hybridization of conduction electrons with 
ocal moments, but as a result of strong antiferromagnetic 
luctuation [l]. Note, that the insulator phase in the state 
)f the Kondo insulator is realized in all compounds with 
)dd number of itinerant electrons per a local moment. 


The Kondo lattice problem is not reduced to a single im¬ 
purity Kondo problem where a local moment is screened at 
low temperatures by conduction electrons and the ground 
state of an impurity is a spin-singlet |4j|5]. Behavior of the 
system (the Kondo problem) is defined by the exponen¬ 
tially small scale (the Kondo temperature) in the exchange 
coupling constant. 

We hope to explain the phenomenon of the Kondo lat¬ 
tice in the framework of the model of a topological Kondo 


lattice proposed. We will call a complex system, which 
includes interacting spin and fermion subsystems one of 
them or two are topological, a topological Kondo lat¬ 
tice. This class of compounds is now known as topologi¬ 
cal Kondo insulators. Recent experiment on SmBg shown 
that a phase state of a topological Kondo insulator is real¬ 
ized in SmB6- The compound SmB6 attracted now atten¬ 
tion due to anomalous behavior of a residual conductivity 
at low temperature T ~ 5 — 7 K that is characterized by 
a 3D to 2D crossover of the transport carriers. The inter- 
metallic compound SuiBg is a strongly correlated insulator 
with topological states at low temperatures, it is the first 
candidate of a topological Kondo insulator family [6-11 . 
According to refs. 12 15 , the topological state in SmBg 
is a result of indirect interaction of 5d-states via strong 
spin-orbit coupling and hybridization between itinerant 
5 d- and narrow (localized) 4/- electronic states. The topo¬ 
logical Kondo insulator is formed via strong interaction of 
dispersive d- and non-dispersive /-electronic bands and 
hybridization between them. Describing a non-dispersive 
band in the framework of the spin operators the system of 
strongly correlated fermions can be defined on the topo¬ 
logical Kondo lattice and the phase of a topological Kondo 
insulator is realized as its the topological state. 

We study a 2D model of the topological Kondo lattice 
at the half-filling and small doping. The model of topolog- 
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Fig. 1: (Color online) The spin sublattice defined on a 
honeycomb-like lattice with links of three directions (x, y 
and z) and z -link broken with x-link; four types of links and 
the re-plaquette operator. 


ical Kondo insulator is the first example of exact solvable 
model of a topological Kondo lattice, the exact solution 
of the Kondo-like lattice model can also be crucial for un¬ 
derstanding of the behavior of the Kondo lattice. In the 
framework of the model we show that the hybridization 
gaps arise in the state of Kondo insulator. The physics of 
the model is surprisingly rich. 

Model. — We focus on 2D model of spin subsystem 
interacting via the Kitaev interaction with a subsystem of 
spinless fermions. The total Hamiltonian 


Hm'Hs+'Hf + 'H, 


int 


(i) 


describes the spin and spinless fermion sublattices, and the 
interaction between them. The Hamiltonian of the spin 
subsystem H s is defined on a honeycomb lattice with a 
spin- \ on each site and an additional exchange interaction 
between spins within the z-links; it can be written in the 
framework of the Kitaev model 
T-links within the z-links 


16 with an additional 


U s = Ax a i a j + Ay °f °j 

( i,j)y 

+ A + ( 2 ) 

where (z, j) is a pair of sites connected by x-, y -, z- and 
T-links mentioned right after the pair notation, aj are 
the three Pauli operators at a site j, and A 7 are the ex¬ 
change integrals along all the links of corresponding di¬ 
rection 7 = x, y, z. The interposed T-link breaks z-link 



Fig. 2: (Color online) The topological Kondo lattice consisting 
of the spin (blue circles) and fermion (red circles) sublattices. 
Red circles form a honeycomb lattice up to a continuous trans¬ 
formation. The interaction between sublattices is defined on 
connected red-blue pairs. 

via the I erf a? exchange interaction (the last term in (l2|) 
with the exchange integral I between spins in two aaai- 
tional sites on the z-link (see Fig. [TJ. For convenience, 
consider A 7 , / > 0. 

For the sublattice of the spinless fermions we will use 
the following tight-binding model that takes into account 
the interaction between nearest and next-nearest neigh¬ 
bors with equal hoppings and pairing amplitudes 

U f = -i*l 52 (d a m + a l a m ) 

( l,m ) 

-it 2 X] ( a \ a rn + a i a mJ + h.C; (3) 

«Z,m» 

where aj and a { are the spinless creation and annihila¬ 
tion operators defined on the honeycomb fermion sublat¬ 
tice and satisfying the usual anticommutation relations, t\ 
is an overlap integral of nearest-neighbor hopping between 
neighbor sites (Z,ra) from l of the type 2 to m of the type 
3, £2 is a clockwise (relatively a honeycomb) next-nearest- 
hopping between second-neighbors ((/,ra)), l and m range 
all sites of both type 2 and 3. 

The unit cell of the lattice consists of four sites of spin 
subsystem (namely 1 — 4 in Fig. [ 2 ]) and two sites of fermion 
subsystem (marked by 2 and 3). The term Hint is gov¬ 
erned by a contact interaction between itinerant states of 
spinless fermions and spin operators at the sites 2 — 3 that 
are interposed into the z-links 

Ui„t = (4) 

l 

where A is the coupling parameter, ni = aja t , l ranges all 
pairs of spin and fermion sites, namely 2 and 3 of each cell 
(see Fig. [2|. 

The Hamiltonian ([!]) defines an exactly solvable model 
of topological Kondo insulator on a honeycomb topological 
Kondo lattice. The Hamiltonian of the spin sublattice H s 
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can be exactly diagonalized using the representation of 
the Pauli operators in terms of a related set of Major ana 
fermions b J and Cj with commutation rules 


{ci,Cj} = 2 Sij, } = 2 5^'Sij, {d, bn = 0, 


with the substitution crj = i tijcj 


16 


(7 = x,y,z). We 
introduce two types of Majorana fermions on each site for 
the sublattice of spinless fermions 


di = a, + a) and gi 



i 



Fig. 3: (Color online) The energy levels of the noninteracting 
(A = 0) fermion subsystem calculated on a cylinder with open 
zig-zag boundary conditions, t\ — 5, £2 = 1 . 


the Hamiltonian PL can thus be rewritten in the form 


U = -1 A ij CiC J ~ dldri 


1/3 

Hj^3 

(3=x,y,z,x (i,j)(3 (Z,m) 


-\t 2 ^2 vimdidm + gidicib \, (5) 


% = "j, 


Dj | 


/phys 


/phys 


with Dj 


bj^jbjCj. 


The operator D 


acts as the identity operator on the physical subspace |16 
it commutes with the Hamiltonian <§• 

The plaquette operator [16] 

defined 


’ sA' 


<%&>° U a2' aV al'°' 


S A" 


IS 


computed the ground-state energy of the model for a set 
of finite size systems and for various set of the exchange 
integrals. In all cases we found that the ground-state en¬ 
ergy is minimized by the same uniform flux pattern. The 
vortex uniform sector with w s = 1 for all plaquette op¬ 
erators w s is the ground state of the model 116, Tt] . The 
model is solved analytically for the uniform configurations, 
due to the translational invariance of the lattice. We now 
focus on the vortex-free configuration (w s = 1 ) for the 
Hamiltonian PL of the entire system ([!]). A contact inter¬ 
action breaks both particle-hole symmetry and time 


A 2 =l, IA = 0, IA = 0 


i 


where the matrix A consists of A]- = A~uJ- for the di- 

l j * i j 

rected links 7 = x,y,z and A x - = Iu x - for the intercalated 


~3 'v 

wise (anticlockwise) next-near-neighbour hopping inside 
corresponding plaquette. 

The physical subspace is defined by the constrain 


IA X =1 
A 2 =0, IA = 0 



A z =0, IA =0 


W S = ^l^S^sA^V 

by a product of the uj- operators around a plaquette s, 

see Fig. |TJ Operators u ^ are constants of motion with 
eigenvalues ±1. Each plaquette operator w s has thus two 
eigenvalues ±1 and it is interpreted as a magnetic flux 
through the plaquette s. The operators ui = i g t b\ are 
the constants of motion with eigenvalues =bl. Thus the 
Hilbert space of states might be split into eigenspaces, 
where all operators u[- and ui might be replaced with 
their eigenvalues. The variables u[- and u\ are identified 
with a static Z 2 gauge fields on the bonds. The Hamil¬ 
tonian © is thus reduced to a quadratic form. The 
itinerant fermions and localized moments reconcile with 
each other via the hybridization with pairing states of 
spin-fermion Majorana fermions on the lattice sites. 

To solve the model exactly, one converts each vortex 
sector to free spinless fermions 16 . Numerically, we have 


Fig. 4: (Color online) Phase diagram of Kitaev’s type for pos¬ 
itive exchange integrals that satisfy IA X + IA y + A 2 Z — 1. 


reversal (TR) symmetry of the model Hamiltonian and 
gives a nontrivial ground-state phase diagram of the sys¬ 
tem. 

Phase analysis. — 

Noninteracting subsystems. In the absence of the in¬ 
teraction @ the system is a sum of two decoupled spin 
and fermion subsystems. Spinless fermion subsystem is 
always in topologically nontrivial phase if ti 7 ^ 0 and 
£2 7 ^ 0 18,19, its spectrum is gapped, one chiral edge 
mode is presented, see Fig. [3] The phase is character¬ 
ized by Chern number C c h of the charge sector, C c h = ±1 


depending on the sign of t\/t 2 - 

Topological phases of noninteracting spin subsystem are 
also associated with Chern number C s . It is convenient to 
illustrate the ground-state phase diagram of the spin sub¬ 
system using Kitaev’s diagram 16 Fig. [4]— a section by 
the plane A 2 Z + IA X + IA y = const in coordinates I A x , 
IA y and A 2 . There are two topological distinct phases: 
topological trivial phase (with gapped states) indexed by 
C s = 0 and topological nontrivial phase (with gapless 
states) with C s = 1 are separated by lines of quantum 
phase transitions. In the case of noninteracting subsys¬ 
tems, an external magnetic field breaks TR symmetry of 
spin subsystem and opens a gap in the spectrum of the 
Majorana fermions [l 6 ] in the gapless region. The struc¬ 
ture of energy levels with edge states in the spin sector is 
shown on Figs [5] for both phases. 

Phase diagram of interacting subsystems. Let us con¬ 
sider the adiabatic connection of the subsystems via the 
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Fig. 5: (Color online) The energy levels of the noninteracting 
(A = 0) spin subsystem calculated on a cylinder with open 
boundary conditions for zig-zag boundary along x-direction, 
gapped 1 = 1 (left) and gapless I — \ (right) case at = 1, 

A x — 2 5 ^y — 2. 


interaction © and evolution of the ground state of the sys¬ 
tem along A and I directions. We use the fixed parameters 
of the Hamiltonian © t\ = 5 and ^2 = 1 in the case of an 
anisotropic spin-exchange interaction \, A y = 2 and 

A z = 1 to demonstrate a band structure and edge modes 
observed in common case. The spectrum of a noninter¬ 
acting fermion subsystem is gapped, a noninteracting spin 
subsystem has gapless spectrum of the Majorana fermions 
when | < \I\ < 1, otherwise it is gapped. There are six 
bands, corresponding to six sites per unit cell: two high 
energy fermion bands and four low energy spin bands (see 
Fig.©. The structure of the edge states depends on a di¬ 
rection of a boundary. We consider the zig-zag boundaries 
along x- and y- exchange bonds for illustration. The direct 
interaction © between subsystems breaks TR symmetry 
of ([!]) and opens the hybridization gap in the spectrum of 
spin-fermion excitations at the half-filling. The interaction 
term © breaks TR symmetry of the complex system, the 
coupling constant A ‘works’ as an external magnetic field 
in the Kitaev model, it opens the gap in the gapless topo¬ 
logical spin phase. The structure of the spectrum and the 
topological number of the entire system do not vary in the 
process of the interaction’s switch and the evolution at the 
adiabatic connection when the value of A increases. The 
Chern number of the system © for arbitrary A depends 
on the Chern numbers of the spin and charge sectors for 
A = 0. There are two possible scenarios of the formation 
of the Kondo insulator phase state along the A direction 
depending of the entry conditions at A = 0: 

• the spin subsystem is in topological trivial phase 
(gapped states), the fermion subsystem is in the state 
of the topological insulator; 

• both the spin and fermion subsystems are in topolog¬ 
ical nontrivial phases. 

Thus the phase diagram of the model system in the co¬ 
ordinates IA x , IA y , IA Z and A is similar to the one in 
Fig.© where the phase diagram of spin subsystem lies in 
the base of a prism and the A axis is normal to the base 
(see Fig. [7]). The planes of the topological phase transi¬ 
tions, in which the bulk gap vanishes in the spin sector 
at A = 0, separate the topological phases with the Chern 
numbers 1 and 2, as it will be shown below. We consider 



Fig. 6: (Color online) The energy levels of the system calcu¬ 
lated on a cylinder with open boundary conditions for zig-zag 
boundary along x-direction with t\ = 5, £2 = 1, A z = 1, 
S= 1, Ay = 2, I = 1 and A = 

( 0 , 0 , 1 ) 



(1,0,0) (0,1,0) 

Fig. 7: (Color online) Phase diagram in coordinates 

(IA X , IA y , A 2 ) with a exchange intergal satisfying a normal¬ 
ization condition IA X + IA y + A 2 = 1 for arbitrary t\, and 
A. 


the cases 1=1 and I = \ to illustrate these two scenarios 
of behavior and investigate peculiarities of the topological 
states. 

Scenario 1. At I = 1 and A / 0 the phase state 
of the system is characterized by two (spin and fermion) 
gaps in the excitation spectrum. Calculations of the Chern 
number and edge states show that C c h = 1 and C s = 0 
at A = 0 , a total Chern number is equal to 1 for arbitrary 
values of A. The structure of edge states depends on a 
direction of folding: ^/-direction boundary (Fig. [ 8 ^l) carries 
chiral gapless edge modes specified by fermion sublattice 
and thus all low-level excitations are fermionic, while x- 
direction boundary (Fig. [ 8 Jd) performs chiral edge modes 
crossing the Fermi level three times with excitations of 
spin type (Fig. [ 8 J 3 ) . These edge modes are associated with 
(7=1, but topological insulator’s behaviors are different: 
it acts as fermion topological insulator in ^/-direction and 
as spin topological insulator in x-direction. 

The interaction term © also directly hybridizes the 
edge states, therefore the edge states are a result of to¬ 
tal spectrum of the system on the one hand, and they are 
a result of direct hybridized interaction on the boundary 
on the other hand. 
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Fig. 8: (Color online) A low energy brushed spectrum (as in 
Fig[6]) at t\ = 5, £2 = 1, I = 1, Aa; = |, A y = 2, A z = 1 with 
boundary along y-direction and a hybridization A = | a) and 
along x-direction and a hybridization A = | b-c). Each point 
has its color depending on fermionic rate of its state which is 
calculated as a projection to the neat fermion state subspace 
(ip\ pr ferm {'ip). A color scale is given. 


Scenario 2. Now consider the second scenario of be¬ 
havior of the Kondo lattice when the closed at A = 0 gap 
in the spectrum of spin excitations opens for A 7 ^ 0 , as 
it is illustrated at I m | in Fig. pi At A = 0 a spin gap 
closing makes the topological index C s ill-defined. Tak¬ 
ing into account a magnetic field as a week disturbance 
that breaks TR symmetry of spin subsystem and opens 
the gap in the spin sector of the spectrum of Major ana 
fermions 16 , we calculate C s and edge states of the spin 
subsystem. According to obtained calculations, C s = 1 
(the same value as C c h)- For A 7 ^ 0, the gap opens due to 
a strong hybridization in a low energy region of the spec¬ 
trum even in the absence of an external magnetic field. 
The total Chern number of the system is equal to two, 
there are two gapless chiral edge modes (hybridized quasi 
spin and fermion ones), see Fig [9| A low energy behavior 
of the edge modes has a spin nature, as we see in Figs |9| 
and as result, a topological spin insulator state with two 
quasi spin edge modes is realized in the Kondo lattice in 
both directions of bounds. 


Overall notes. The model of the proposed 2 D topo¬ 
logical Kondo lattice takes into account an interaction be¬ 
tween different topological subsystems, namely, spin and 
fermion lattices. The spin-fermion interaction leads to the 
hybridization of the Majorana fermions of the subsystems. 
The considered interaction is reduced to the hybridiza¬ 
tion between excitations in the spin and charge sectors, 
as result the Chern number of the system is characterized 
by one (total) Chern number, that depends on spin and 
charge Chern numbers of the corresponding noninteract¬ 
ing subsystems. 

In spite of the Hamiltonian ([5| does not conserve the 
total number of particles, the Fermi surface is sharply 
defined. The hybridization of itinerant spinless fermions 


Fig. 9: (Color online) A low energy brushed spectrum (as in 
Fig[8]) at t± = 5, t 2 = 1, / = §, A = = §, A y = 2, 

A z = 1 with boundary along x-direction a) and y-direction 
b-c). A color scale is given. 


and local moments form the volume of the Fermi sea 
f BZ dk = n + 1 , where n is a density of fermion and 
integrating is over the Brillouin zone of the system. A 
stability of the phase of the Kondo insulator does not de¬ 
pend on the eigenvalues of the local operators ui = 2 tp — 1 
because the ground state energy is realized both tp = 1 
and rii = 0 , where tp is a number of fermions on the site 
l of the fermion sublattice. Therefore there are free bulk 
states of Majorana fermions (with zero energy) located at 
the sites of the fermion sublattice. The total ground-state 
energy of the system reaches minimum at the half-filling 
of fermion and spin subsystems n — 1 , this phase state 
is also realized under an arbitrary doping p = n — 1 of 
particles, because the rest of p particles are located at the 
fermion sublattice (on the sites with tp = 1 ) and occupy 
the Fermi surface ef = 0 in the insulator phase. The 
total energy of the system remains the same, it does not 
depend on the doping of the system, a static Z 2 gauge field 
‘works’ as a reservoir for doped particles. In contrast to 
traditional insulator phase, the Fermi surface is occupied 
and an insulator gap is effectively two times less. 

Conclusions. — In summary, we have considered the 
implication of the adiabatic connection between spin and 
fermion subsystem defined on a honeycomb Kondo lattice. 
There exist two possibilities of realization of the Kondo in¬ 
sulator state, both cases have been considered. We have 
calculated a hybridization gap in the framework of the 
model proposed (up to now the mechanism of forming of 
the hybridization gap in the state of the Kondo insulator 
has been unknown). The spinless fermions are localized 
at the lattice sites due to the contact interaction with lo¬ 
cal moments, this phenomenon is analogous to the Kondo 
screening in real space (on lattice sites). The itinerant 
states of spinless fermions hybridize with spin excitations, 
as a result, the gap opens at the half-filling. The gapless 
edge states form a surface subband of chiral Majorana 
fermions. It is shown that hybridization within topologi- 
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cal Kondo insulator can lead to changing fermionic topo¬ 
logical insulator into spin one in sense of low energy edge 
excitations. 

In the spin-wave approach, the spin excitations are 
bosons, they are not hybridized with fermions. The 
ground state of the Kondo lattice at half-filling can be 
explained by a hybridization gap that arises from a hy¬ 
bridization of itinerant electrons and local moments. In 
this letter we have firstly proposed the mechanism of 
hybridization between local spins and itinerant spinless 
fermions in the framework of the proposed exactly solv¬ 
able model. 
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